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There are two physically different interpretations of "triviality" in (A$^)4 theories. The 
conventional description predicts a second-order phase transition and that the Higgs mass 
rrih must vanish in the continuum limit if v, the physical v.e.v, is held fixed. An alternative 
interpretation, based on the effective potential obtained in "triviality-compatible" approx- 
imations (in which the shifted 'Higgs' field h{x) = <l>(x) — ($) is governed by an effective 
quadratic Hamiltonian) predicts a phase transition that is very weakly first-order and that 
rrih and v are both finite, cutoff-independent quantities. To test these two alternatives, we 
have numerically computed the effective potential on the lattice. Three different methods 
were used to determine the critical bare mass for the chosen bare coupling value. All give 
excellent agreement with the literature value. Two different methods for obtaining the ef- 
fective potential were used, as a control on the results. Our lattice data are fitted very well 
by the predictions of the unconventional picture, but poorly by the conventional picture. 



1 Introduction 



One of the most interesting results of modern quantum field theory concerns the "triviality" 
Jl], 1^, 1^, §, H, H, of (A$^)4 theories. The physical meaning of this mathematical result re- 
mains controversial, however. The conventional interpretation is based on Renormalization- 
Group- Improved-Perturbation-Theory (RGIPT), while a quite different interpretation is ad- 
vocated in Refs. ||, ^, |Ty]. The two pictures have drastically different imphcations for the 
Standard-Model phenomenology. They also give different predictions for the effective poten- 
tial, and the purpose of this paper is to compare those predictions with a model-independent 
lattice calculation of Vcs- That is, we perform a precise numerical experiment as a test of 
the two alternatives. 

The conventional interpretation of "triviality" is as follows: Leading-order RGIPT pre- 
dicts that the running coupling constant, if finite at some low energy scale, will blow up 
and become infinite at some larger energy, the Landau pole. The only way to avoid this 
unphysical behaviour, it is argued, is to push the Landau scale to infinity by sending the 
low-energy renormalized coupling A^j to zero, thus recovering "triviality." In RGIPT the 
Higgs mass rrih is proportional to (Xrv'^Y^'^ and goes to zero in the continuum limit, if the 
vacuum expectation value v (phenomenologically determined to be ~ 246 GeV) is taken to 
be finite. In this picture, the only way to have a viable Higgs mass is to keep the cutoff A fi- 
nite and not too large. The scalar sector of the Standard Model is then a non-renormalizable 
theory in which the mysterious cutoff plays a crucial role. 



One problem with this explanation is that RGIPT does not give a consistent picture |TT 
The Landau pole appearing in leading order is absent in next-to-leading order. Instead, 
there is an ultraviolet fixed point A, due to the negative sign of the 2-loop coefficient of 
the perturbative /3-function. Taken at face value, the two-loop result implies a finite bare 
coupling constant, with the renormalized coupling lying anywhere in the region to A. Since 
there is no reason for A/j to vanish in this case, one has a direct conflict with "triviality" 
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It is usually asserted that the leading-order RGIPT picture is supported by non-perturbative 
lattice simulations showing that ruh vanishes in units of v in the continuum limit. While the 
evidence certainly implies rrih/vB — »• (in accord with the rigorous results of sect. 15 of 0] 
that rrih and Vb cannot scale uniformly in the continuum limit), the crucial issue is how the 
bare vacuum field vb measured on the lattice is related to the physical v ~ 246 GeV defined 



from the Fermi constant. Up to now [|I^ one has used the long-distance behaviour of the 
propagator of the shifted field(s) to extract a renormalization constant Z = Zh, and defined 
V = VB/\/Zh- The lattice data provide overwhelming evidence that Z^ ~ 1, as one would 
expect for a trivially free shifted field. However H, |], |lOl, is VB/\/Zh the correct definition 
of the physical vacuum field? More precisely, if one considers the exact definition of the 



physical vacuum field from the effective potential, namely (pR = 0^/^/^^ such that 



^1 



ml, 



--±v 



does one find Z^ = Zh up to negligible corrections? 

In fact, as discussed in refs. [|], ^ [1^], one finds Z^f) ^ Zh in any approximation (e.g., 
one-loop, gaussian p, or post-gaussian [|1^) that mimics the basic "triviality" of the the- 
ory, where the shifted field h{x) = $(x) — 0^ is consistently described by an effective 
quadratic Hamiltonian. In such approximations there is a non-trivial Z^, although Zh = 1 
holds identically. The resulting effective potential, Kriv, has a simple, universal structure 
given by the sum of a (suitably redefined) classical potential and a (suitably redefined) 
zero-point energy for a free field with a 0-dependent mass. ["Suitably redefined" refers to 
the usual mass-renormalization or normal-ordering procedure.] The crucial point is that 
this simple structure - originally obtained in the Coleman- Weinberg one-loop calculation 
15[ - is a nonperturbative consequence of "triviality". The fact that the one-loop approx- 



imation seems untrustworthy from a loop-expansion perspective |T^ is not relevant. The 
one-loop computation should rather be viewed as the prototype for a class of nonperturba- 
tive, "triviality-compatible" computations of Ves, which all yield the same result, Kriv The 
form of Vtriv reflects the coexistence of spontaneous symmetry breaking (SSB) and "trivial- 
ity", and becomes the basis for an alternative renormalization procedure that replaces the 
standard perturbative approach. 

To understand the basic difference with the usual perturbative approach, consider the 
following question: In the continuum limit of a "trivial" theory, where the S-matrix for 2 — i> 2 
particle scattering reduces to the trivial identity, can one still generate a finite energy density 
to de-stabilize the perturbative vacuum? Indeed, yes, and statistical physics is full of such 
examples; e.g. superconductivity, where an arbitrarily small 2-body interaction produces 
macroscopic effects. The physical mechanism is well known: a tiny 2-body interaction g can 
produce drastic changes in the vacuum structure if there is a sufficiently large number of 
states at the Fermi surface. In this situation, where the energy gap in the single-particle 
spectrum is a many-body effect, ordinary perturbation theory in g fails in predicting all the 
basic features of the superconducting ground state. 

Following this line of thought one deduces the simple physical picture of refs. p, |^, p!0|] 
where the (nearly) massless quanta of the symmetric phase condense in the zero 4-momentum 
state. Even with an infinitesimal 2-body strength, the Bose condensation produces a finite 
gain in the energy density, leading to the instability of the perturbative vacuum. The 
excitations of the new vacuum are non-trivially related to the original quanta, but they 
also have vanishingly small interactions. Therefore, a proper renormalization procedure for 
(A$^)4 theory cannot be based on a vain attempt to generate a finite 'Ar', a concept which 
has no place in a "trivial" theory. Rather, the correct strategy has to be based on the 



physical requirement that the energy density associated with SSB is finite [|], |J, |T6|, |1^ |1S 



191 , |20| , pT| , p^ ; i.e. one should require the effective potential to be cutoff independent. A 
straightforward renormalization-group analysis of Vtnv then shows that, while miJvB —^ 0, 
one gets irth/v = cutoff-independent. Clearly, this has radically different implications for 
the Standard Model. The non-trivial re-scaling of the vacuum field, v = vb/ ^/Z^, where 
Z^ — >■ oo, is quite distinct from the trivial Z^ = 1 renormalization of the free fluctuation 
field. This is the essential ingredient that represents in a quantum-field-theoretical context 
the intuitive notion of an infinitely dense ||2^ Bose condensate at p^ = coexisting with 
trivially free excitations at p^ ^ 0. 



Another striking feature of V^riv is that it predicts a first-order phase transition p|, |T^. 



By contrast, the RGIPT result, V^crt, shows a second-order transition. With Vtriv one sees 
that, as the bare-mass term tq = m?^ is made more and more negative, the SSB transition 
occurs at a value tq = r^ where m, the physical mass-gap of the symmetric phase, though 
infinitesimal in units of the corresponding mh, is still non-zero. Only at an even more 
negative value of the bare mass, tq = r^ < r^, does one find the 'Coleman- Weinberg regime' 
where m = identically. However, by then the system is well inside the broken phase. Since 
^c 7^ "^s there is no continuum limit for non-zero values of the bare coupling Aq, contrary 
to the prediction of leading-order RGIPT where m always vanishes exactly at the phase 
transition. Only if Aq — > can one obtain a continuum limit from Kriv One finds that 
\{tc — ''"s)/'"c| ~ exp(— 87r^/3Ao) 0, so that r^ — > r^ as Aq — > 0, yielding a transition that 
asymptotically becomes second order, in agreement with the rigorous result known for the 
gaussian model. The difference between r^. and r^, although exponentially small for weak 
bare coupling, remains crucial |1^. Indeed, in the limit, infinitesimal vanations of tq near 



the phase-transition value induce finite variations in the particle mass of the broken vacuum; 
in the energy-density difference between the two phases; and in the barrier between the two 
phases. The problem with the conventional approach is that it looks at the phase transition 
on too coarse a scale - making finite variations in tq. Viewed on that scale the transition 
appears indistinguishable from a second-order phase transition and the fine details are not 
seen. 

The qualitative difference between Vi;riv and l^ert means that the two pictures can be 
distinguished by a sufficiently precise lattice calculation of Ves- That is, one can perform a 
model-independent, numerical experiment to test the predictions of both the conventional 
RGIPT picture and the alternative picture of refs. ||, |], |10| . Initial results were presented 



in ref. |2^ . The aim of this paper is to provide more refined results from a precise lattice 



calculation of the slope of the effective potential. 



2 The lattice effective potential 



We begin by defining the (one-component) (A$^)4 theory on a lattice: 



S = a*Y, 






where x stands for a generic lattice site, a denotes the lattice spacing, and Aq > 0. For SSB 
the basic quantity is the expectation value of the bare scalar field $(2;) (B=Bare) 

($)j = MJ) (2) 

in the presence of an external source whose strength J is x-independent. Determining 0_b( J) 
at several J-values is equivalent [0, |2^ to inverting the relation 



J = J{<f>B) = ^ (3) 

d(pB 

involving the effective potential Ves{4'B)- In this way, starting from the action in Eq. (|l|), the 
effective potential of the theory is rigorously defined up to an arbitrary integration constant 
(usually chosen to fix KflF(O) = 0). [This definition is equivalent to the Legendre transform 
definition and is convex downward p^.] In this framework, SSB occurs when the function 
4>b{J) has a non-zero limit as J ^ 0: 

jnn 0b(J) = ±VB ^ 0. (4) 

One expects such behaviour for a certain range of the bare parameters tq and Aq appearing 
in the lattice action Eq. (|l]). It corresponds to the effective potential having non-trivial 
minima with Vcs{^vb) < V'cff(O). 



2.1 Monte Carlo simulation 

For the Monte Carlo simulation of the lattice field theory described by Eq. (|I]) we used the 
standard Metropolis algorithm. In order to avoid the trapping into metastable states due 
to the underlying Ising dynamics we followed the upgrade of the scalar field $(a;) with the 
upgrade of the sign of $(a;). This is done according to the effective Ising action [p7| 



5'ising = J5^|$(a;)|s(x) -^^|<l>(x + /i)$(x)|s(x + /i)s(x) , (5) 

X X fl 

where s{x) = sign($(a;)). We measured the vacuum expectation value of the scalar field 

1 ^^ 1 



1=1 



where N^. is the number of the lattice configurations generated with the action, Eq. (|l]). 

Statistical errors are evaluated taking into account the autocorrelation time in the sta- 
tistical sample generated by the Monte Carlo simulation. If we consider a generic observable 
O (function of the lattice configuration) the integrated autocorelation time is defined as 

_ 00 _. 00 

rUO) = 2 E P<^o{t) = 2 + E Pooit) , (7) 

t = —00 t = 1 



where poo{t) is the normahzed autocorrelation function: 

Poo{t) = Coo(t)/Coo(0), (8) 

and 

Coo{t) = {OsOs+t) - {Of (9) 

is the unnormahzed autocorrelation function. 

The integrated autocorrelation time depends on the parameters in the lattice action 
Eq. (|I|) and determines the statistical error in Monte Carlo measurements of the expectation 
value {O) of the observable O. If O is the sample mean of the observable O 

1 " 

= - 5^ o*, (10) 

the sample variance is 



:{0) = m - {Of = \J2 Cooir - .) = - J2 (^--) 

r,s=l t = -(n-l] ^ ^ 



^ - (2r,„i(0)) Coo(O) for n > r . (11) 

n 



Therefore the statistical error is given by 



s=\\-Coo{^). (12) 

V n 

For completely uncorrected data K = 2Tint{0) = 1. 

Therefore, to obtain the statistical error we must estimate the factor K in Eq. (|12|). This 



can be achieved through a direct evaluation of the integrated autocorrelation time pSf , or by 



using the "blocking" |^ or the "grouped jackknife" [^] algorithms. We have checked that 
applying these three different methods we get consistent estimates of the statistical errors. 



2.2 Determination of the critical bare mass parameter 

We chose to run our lattice simulation with Aq = 0.5. We then have s = 3Ao/167r^ -C 1, as 
needed for the continuum limit of refs. |§|, ^, . This puts us in a region where both bare 
and renormalized couplings are small: an excellent place to test the validity of perturbation 
theory. 



The bare mass-squared tq needs to be close to the critical value Vc so that the correlation 
length will be very large. To determine accurately the r^ for Aq = 0.5 we used several 
methods. Firstly, there is an analysis by Brahm |^l| that yields 



1.00 ±0.05 , , 

re = -0.2240 ±0.0010, (13) 

for an L^ size lattice. For L=16 this gives 

re = -0.2279(10) . (14) 

Brahm made use of the Liischer-Weisz high-temperature-expansion results |^ and made 
lattice calculations of the susceptibility in the broken phase on lattices ranging from 4*^ to 
8^. The susceptibility x is defined as: 

where 

X 

One expects that near the critical region x^^ ~ (r^ — Tq), modulo logarithmic corrections 
to the free-field scaling law. One can thus determine r^ by extrapolation to vanishing x~^- 
Strictly speaking, this method is valid only for a second-order phase transition where r^ = r^ 
and both m and rrih vanish at the phase transition. In the case of a very weak first-order phase 
transition where \rc — r^l/rc ~ exp(— l/(2s)) the induced uncertainty should be negligible. 

To check Brahm's result we performed three different numerical calculations of Tc at 
Ao = 0.5 on a IG'^ lattice. We determined Tc from the susceptibility Eq. (p!5| ) in both the 
broken and symmetric phases. Our data are well described by the simple linear fit 

X~'^ = a\r - rd . (17) 

We find in the broken phase: 

a = 1.964(40) 
r^ = -0.2270(14) 
xVd.o.f. = 0.79, (18) 

and in the symmetric phase 

a = 2.529(70) 
re = -0.2296(16) 
xVd-o.f. = 0.33. (19) 



Our data do not show evidence of the logarithmic corrections. As a matter of fact we also 
tried the fits: 

X^'^ = a\r — rd I In \r — VcW^ , (20) 

and found 7 consistent with zero in both cases [0.020±0.121 (broken phase), — 0.003±0.133 
(symmetric phase)]. Our results for x~^ together with the fits Eqs. (0), (|1^) are displayed 
in Fig. 1. 

Our third calculation of Vc was obtained from the generalized magnetization ($) which 
should have the form [R2i: 



($) = a{rc - ry^^\\n\r - r^\f + 6 . (21) 

Accordingly we fit our data and found 

a = 1.459(17) 

re = -0.2278(19) 

p = 0.0064(187) 

S = 0.0095(169) 

xVd-o.f. = 0.27. (22) 



The data and the fit are shown in Fig. 2. Combining the various estimates in Eqs.(|T8|, |T£ 
22) our final value of Vc is: 



r^ = -0.2280(9) . (23) 

The agreement with Eq. (|T4[) is excellent. We thus have three independent confirmations 
that the Brahm's value Vc = —0.2279, extrapolated from smaller lattices by means of Eq. 
([T3|), represents a precise input definition of the 'Coleman- Weinberg regime' on a 16^ lattice 
with the action Eq. (1) at Aq = 0.5. 



2.3 Determination of the effective potential 

We have used two independent methods to compute the effective potential. Firstly, we ran 
simulations of the lattice action Eq. (^ for 16 different values of the external source in the 
range 0.01 < | J| < 0.70. In this way, as outlined in Eqs. (2-4), we directly obtain the slope 
of the effective potential (from which Ves can be obtained, up to an additive integration 
constant). Our results for ($)j = (psiJ) are shown in Table 1 (errors are statistical only). 

As an additional check of our results, we performed a calculation using an alternative 



approach to V^s first proposed in ref. [^ and later extended in ref. |^. This second 



J 


MJ) 


J 


MJ) 


-0.010 


-0.288862 (695) 


0.010 


0.289389 (787) 


-0.030 


-0.413565 (321) 


0.030 


0.414713 (376) 


-0.050 


-0.488797 (296) 


0.050 


0.489132 (249) 


-0.075 


-0.557737 (181) 


0.075 


0.557961 (182) 


-0.100 


-0.612497 (169) 


0.100 


0.612865 (151) 


-0.300 


-0.876352 (111) 


0.300 


0.876518 (95) 


-0.500 


-1.03526 (8) 


0.500 


1.03532 (7) 


-0.700 


-1.15518 (8) 


0.700 


1.15528 (7) 



Table 1: Values of 0_b(J) obtained on a 16"^ lattice at Aq = 0.5 and tq 
Errors are statistical only. 



-0.2279. 



method is based on the approximate effective potential Uesi<pB', ^) defined through 

exp -f/eff(0B; Q) = [[D^ S (^ fd^x <l>(x) - ^b) exp -S[$]. (24) 

In the limit in which the 4-volume fi — ;> oo, Ues tends to the exact Kff(0s) from Eqs. (||, 
H). The difference between Ues{(pB',^) and Ves{(pB) gives both a consistency check of our 
calculations and an indication of the effects due to the finiteness of our lattice. For our 
purposes, it is more convenient to compare 7(0^), from Table 1, with the corresponding 
quantity in the alternative method p4[ 



Jeff(0B,fi) 



dUcs{(pB;^) 



■>B 



^o{'^^)^B +ro(f)L 



where 



{^%,^^{Jd'x^'{x)) 



(25) 



(26) 



and all expectation values (...) are computed holding ($) = (pB = fixed. 

Let us give some more details about this second approach. The constraint on the value 
of ($) is implemented by updating a pair of sites at the same time so that ($) = 0^ remains 
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0B 


Jes{4>B; fi) 


(pB 


Jcs{(J)b;^) 


-0.288862 


-0.00976 (7) 


0.289389 


0.00986 (9) 


-0.413565 


-0.02983 (7) 


0.414713 


0.03008 (9) 


-0.488797 


-0.04992 (10) 


0.489132 


0.04999 (8) 


-0.557737 


-0.07485 (8) 


0.557961 


0.07493 (8) 


-0.612497 


-0.09980 (9) 


0.612865 


0.09995 (8) 


-0.876352 


-0.29978 (12) 


0.876518 


0.29992 (10) 


-1.035260 


-0.49991 (12) 


1.035320 


0.49993 (11) 


-1.155180 


-0.69970 (15) 


1.155280 


0.69992 (14) 



Table 2: Values of Je{i{(pB] ^) from Eq. (p5D, as obtained with our 16^ lattice at Aq 
and ro = r^ = -0.2279. 



0.5 



constant. Then we compute the action variation for this double change and accept or reject 
it by Metropolis algorithm. The generation of the sites pair is such that at least one member 
of the pair sequentially spans the whole lattice. Afterwards, we perform an Ising update of 
the field signs (by using Eq. (|^) |^ ) to avoid unwanted trapping also in this case. After a 
run we use the jackknife algorithm ^^ to evaluate {^^)^g with its statistical error and use 
Eq. (^) to find the value of Jes and its associated statistical error. 

The results from this computation of Jes{(pB',^) are reported in Table 2. The input 
values of 0^ were chosen to be the output values from Table 1. 

To get the total statistical errors reported in Table 2 we have combined in quadrature the 
purely statistical error on JcAF at any fixed value of 0^ with that obtained by propagating the 
errors on (pB reported in Table 1. To this end we have used a fitting function which provides 
an excellent fit to the data J = J{(f)B) from Table 1. This estimate of the total statistical 
error has been checked for consistency by performing, for a few values of 0b, two runs at 
(j)B ± S(f)B, S<pB being the statistical error affecting 0^ as determined from Table 1. The two 
estimates give essentially equivalent results. 



3 Comparing theory with the lattice data 

We can now compare the lattice data in Tables 1 and 2 with the existing theoretical expecta- 
tions. In the case of refs. [|], [10[, the predicted form (in the Coleman- Weinberg case, tq = Tc, 
where no quadratic term is present in the effective potential) is: 



^tr\\ 



dVtri^ 



'B 



a0|ln(0|)+70|, 



(27) 



JB 



where a and 7 are free parameters. (Their values are approximation-dependent within the 
class of "triviality-compatible" approximations.) 

The RGIPT prediction exists in various slightly different forms in the literature. We have 



first used the full two-loop calculation of Ford and Jones |^ in the dimensional regularization 
scheme. Their expression for the effective potential, in a one-component theory, is (A = 6A0) 



V 



2— loop / 



with 



^B 



Vr 



oyvB, 



+ V, 



IW^B, 



+ V2' 



>B 



{2i 



VMb) 



^4 + —01 



1 4^m^ 3 

-m<^\ul — 

647r2 2L ^2 2^ 



and 



V, 



2WB 



X^(f)%mi 



2567r4 8 



-— mn 



5 + 8^(1) - 41n^ + In 



:r-2Tnn 



fl^ 



+ 



Xm2 



2567r4 8 



/i2 



-\ 2 



• (29) 



In the above equations we have introduced m^ = A0^/2 + M^, fl{l) = jS — |C(2) with 
S = 1/2^ + 1/5^ + 1/8^ + . . . while In includes in the definition of the logarithm additional 
terms of the MS scheme. The prediction for J follows by differentiation: 



t2— loop/ 



^^2-loop 



^B 



'JB 



(30) 



In this case, the two free parameters are the scale /x, and the mass parameter M^ of the 
classical potential. 

A different version of the RGIPT prediction, which re-sums various terms, is given in 
Eq. (242), Sect. 5.4.2, of the textbook by Itzykson and Drouffe (ID) P§, namely 



^'^(0b) 



Act) 



B 



In J^ 



(4vr)V| 



V3 i81nJf^ 
\<pb\ 



+ 



(31) 
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data 


J*"" Eq. (p7|) 


J2-100P Eq_ (|3Qj) 


JID Eq. (31) 


Table 1 


a = 0.0152(2) 
7 = 0.4496(1) 
y2 - 15 

'^ 16-2 


/i = 8.0304(449) 
M2 = -0.0025(1) 

2 _ 142 
^ 16-2 


|;U = 2.70(1) X 10^ 
A = -0.0055(3) 

2 _ 116 
■^ 16-2 


Table 2 


a = 0.0156(2) 
7 = 0.4494(1) 

2 _ 10 
'^ 16-2 


/i = 7.9883(455) 
M2 = -0.0028(1) 

2 _ 109 
'^ 16-2 


Ifi =2.69(1) X 10^ 
A = -0.0063(3) 

2 _ 85 
^ 16-2 



Table 3: We report the values of the parameters together with the x^ obtained by fitting 
Eqs. (|27|) , (|30|) , and (|3lD to the data reported in Table 1 and 2. 



where again we have two free parameters A and /i. 

Note that in the two preceding equations we have ignored the distinction between (pB 
and (pFt- This is justified because in these conventional approaches there is only one Z (i.e., 
Z^ = Zh) and it is known from many lattice calculations that Zh is very close to unity (see 
TableIIofref.[|13l). 



The results of fitting the lattice data to the three predictions, Eqs. ( P7| , PO] , PT| ) are 
reported in Table 3. The first yields a good fit (x^/d.o.f. < 1), while the latter two yield poor 
fits (x^/d.o.f. ~ 6 — 10). Thus, the data significantly favour the unconventional interpretation 
of "triviality" proposed in refs. H, |^, [1^] over the conventional interpretation. 



4 Conclusions 



In this paper we have presented a numerical experiment to test the two alternative pictures 
of 'triviality' presented in the Introduction. To this end we have first determined the value 
of the critical bare mass parameter r^ = r^ that defines the 'Coleman- Weinberg- regime' 
of (A$*i)4 theory on our 16^ lattice for Aq = 0.5. Using three different methods we have 
confirmed the pre-existing estimate r^ = —0.2279 obtained by Brahm |31]. Then, we have 
computed the effective potential, using two different methods as a control on our results. 
The quality of the fits to the lattice data is important evidence for the unconventional 
interpretation proposed in [|8|, |9|, |10 . 
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FIGURE CAPTIONS 

Figure 1. The inverse of the susceptibihty Eq.(^) vs. the bare mass-squared Vq, fitted by Eq.(|T 
in the symmetric (left side) and broken-symmetry phase (right side). 

Figure 2. The generahzed magnetization $ vs. the bare mass-squared Tq, fitted by Eq.(pTD in 
the broken-symmetry phase. 
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Figure 1: The inverse of the susceptibihty Eq. ([l5|) vs. the bare mass-squared vq, fitted by 
Eq.(|l^) in the symmetric (left side) and broken-symmetry phase (right side). 
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Figure 2: The generalized magnetization $ vs. the bare mass-squared tq, fitted by Eq.(|2l] 
in the broken-symmetry phase. 
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